Coherent Backscattering of Ultra-cold Atoms 
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We report on the direct observation of coherent backscattering (CBS) of ultra-cold atoms, in a 
quasi two dimensional configuration. Launching atoms with a well defined momentum in a laser 
speckle disordered potential, we follow the progressive build up of the momentum scattering pattern, 
consisting of a ring associated with multiple elastic scattering, and the CBS peak in the backward 
direction. Monitoring the depletion of the initial momentum component and the formation of the 
angular ring profile allow us to determine microscopic transport quantities. The time resolved 
evolution of the CBS peak is studied and is found a fair agreement with predictions, at long times 
as well as at short times. The observation of CBS can be considered a direct signature of coherence 
in quantum transport of particles in disordered media. It is responsible for the so called weak 
localization phenomenon, which is the precursor of Anderson localization. 

PACS numbers: 03.75.-b, 67.85.-d, 05.60.Gg, 42.25.Dd, 72.15.Rn 



Quantum transport differs from classical transport by 
the crucial role of coherence effects. In the case of trans- 
port in a disordered medium, it can lead to the com- 
plete cancelling of transport when the disorder is strong 
enough: this is the celebrated Anderson localization 
(AL) pp. For weak disorder, a first order manifestation of 
coherence is the phenomenon of coherent backscattering 
(CBS), i.e. the enhancement of the scattering probability 
in the backward direction, due to a quantum interference 
of amplitudes associated with two opposite multiple scat- 
tering paths [2]-[4] (see inset of Fig. [I]). Direct observation 
of such a peak is a smoking gun of the role of quantum 
coherence in quantum transport in disordered media. 

CBS has been observed with classical waves in op- 
tics [5 - 8 , acoustics [9j [TOj , and even seismology [11] . In 
Condensed Matter Physics, CBS is the basis of the weak 
localization phenomenon (see e.g. [12 j), which is respon- 
sible for the anomalous resistance of thin metallic films 
and its variation with an applied magnetic field [13]. In 
recent years, it has been possible to directly observe An- 
derson localization with ultra-cold atoms in one dimen- 
sion p~4] [15] and three dimensions [16] [17] . Convincing 
as they are, none of these experiments includes a direct 
evidence of the role of coherence. 

In this letter, we report on the direct observation of 
CBS with ultra-cold atoms, in a quasi two dimensional 
(2-d) configuration [15] . Our scheme is based on the pro- 
posal of Ref. [20 that suggested observing CBS in the 
momentum space. A cloud of non-interacting ultra-cold 
atoms is launched with a narrow velocity distribution in a 
laser speckle disordered potential (Fig.[T]). Time of flight 
imaging, after propagation time t in the disorder, directly 
yields the momentum distribution, as shown in figure [2] 
As expected for elastic scattering of particles, we observe 
a ring that corresponds to a redistribution of the momen- 




FIG. 1. Experimental set-up. A cloud of non-interacting 
ultra-cold atoms, released from an optical trap (red beams) 
and suspended against gravity by magnetic levitation (yellow 
coils), is launched with a well defined momentum pi along 
the z-axis. It is submitted to an anisotropic laser speckle dis- 
ordered potential (blue beam), elongated along the x axis, 
leading to a quasi 2-d diffusive motion in the y — z plane 
(see text). The atomic momentum distribution is monitored 
by fluorescence imaging after a the time of flight of 150 ms. 
Inset: physical origin of CBS. The coherent enhancement of 
scattering in the backward direction originates from the in- 
terferences between each multiple scattering path (solid line) 
and its reversed counterpart (dashed line). 



turn directions over 2tt while the momentum magnitude 
remains almost constant. The evolution of the initial 
momentum peak and of the angular ring profile yields 
the elastic scattering time and the transport time. But 
the most remarkable feature is the large visibility peak, 
which builds up in the backward direction. The height 
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and width of that peak, and their evolution with time, 
are an indisputable signature of CBS, intimately linked 
to the role of coherence. 

To understand the origin of that CBS peak, let us con- 
sider an input plane matter-wave with initial momentum 
Pi that experiences multiple scattering towards a final 
momentum pf (inset of figure [T]). For each multiple scat- 
tering path, we can consider the reversed path with the 
same input pi and output pf. Since the initial and final 
atomic states are the same, we must add the two cor- 
responding complex quantum amplitudes, whose phase 
difference is 5<j) = (pi + pf) • R/ft (R is the spatial sepa- 
ration between the initial and final scattering events and 
h = h/2ir the reduced Planck constant). For the ex- 
act backward momentum pf = — pi, the interference is 
always perfectly constructive, whatever the considered 
multiple scattering path. This coherent effect survives 
the ensemble averaging over the disorder, so that the to- 
tal scattering probability is twice as large as it would 
be in the incoherent case. For an increasing difference 
between pf and — pi, the interferences are progressively 
washed out as we sum over all interference patterns asso- 
ciated with all multiple scattering paths. It results in a 
CBS peak of width inversely proportional to the spread 
AR in the separations [21]. For diffusive scattering 
paths, the distribution of R is a Gaussian and the CBS 
widths decrease in time according to ApcBS,a = h/y/2Dt 
for each direction of space (a = y,z), D being the dif- 
fusion constant along that direction. This time resolved 
dynamics of the CBS peak has been observed in acous- 
tics [11[10J and optics |22ll23j. 

The crux of the experiment is a sample of non- 
interacting paramagnetic atoms, suspended against grav- 
ity by a magnetic gradient (as in [H]), and launched 
along the z-axis with a very well defined initial mo- 
mentum pi (see Fig. [TJ. This is realized in four steps. 
First, evaporative cooling of an atomic cloud of 87 Rb 
atoms in a quasi-isotropic optical dipole trap (trapping 
frequency ~ 5 Hz) yields a Bose Einstein condensate 
of 9 x 10 4 atoms in the F = 2, mp — — 2 ground 
sublevel. Second, we suppress the inter-atomic interac- 
tions by releasing the atomic cloud and letting it ex- 
pand during 50 ms. At this stage, the atomic cloud 
has a size (standard half-width along each direction) 
of Ar a = 30 /im, and the residual interaction energy 
(E[ nt /h ~ 1 Hz) is negligible compared to all relevant 
energies of the problem. Since the atomic cloud is ex- 
panding radially with velocities proportional to the dis- 
tance from the origin, we can use the "delta-kick cooling" 
technique [25 , by switching on a harmonic potential for 
a well chosen amount of time. This almost freezes the 
motion of the atoms, and the resulting velocity spread 
Av a = 0.12 ± 0.03 mm/s is just one magnitude above 
the Heisenberg limit (Ar a • mAv a ~ bh). Last, we give 
the atoms a finite momentum pi along the z-direction, 
without changing the momentum spread, by applying an 



A 0.4 




Pz 




^ ms 




FIG. 2. Observed momentum distributions after different 
propagation times t in the disorder. The images correspond to 
an averaging over 20 experimental runs. Note that the scale 
is different in the three first images (t =0, 0.5, and 1 ms), 
whereas it is the same in the three last images (t =1.5, 2, and 
2.5 ms). 



additional magnetic gradient during 12 ms. The first 
image of Fig. [2] shows the resulting 2-d momentum dis- 
tribution. The average velocity is v x = 3.3 ± 0.2 mm/s 
(k[ = Pi/h — 4.5 /im -1 ), corresponding to a kinetic en- 
ergy E K = pf /2m (E K /h ~ 1190 Hz). This momentum 
distribution is obtained with a standard time of flight 
technique that converts the velocity distribution into a 
position distribution. Because of the magnetic levita- 
tion, we can let the atomic cloud expand ballistically for 
as long as 150 ms before performing fluorescence imag- 
ing along the x axis. The overall velocity resolution 
Av res = [Av a 2 + (Ar a /t to f) 2 ] 1/2 = 0.23 mm/s, is nev- 
ertheless mainly limited by the size Ar a of the atomic 
cloud. 

To study CBS, we suddenly switch on an optical dis- 
ordered potential in less than 0.1 ms, let the atoms scat- 
ter for a time t, then switch off the disorder and mon- 
itor the momentum distribution at time t. The disor- 
dered potential is the dipole potential associated with 
a laser speckle field [26j [27] , obtained by passing a laser 
beam through a rough plate, and focusing it on the atoms 
(Fig. [I}. It has an average value Vr (the disorder "am- 
plitude") equal to its standard deviation. Its autocor- 
relation function is anisotropic, with a transverse shape 
well represented by a Gaussian of standard half-widths 
&y — o z = a± c± 0.2 /im, and a longitudinal Lorentzian 
profile of half- width a x ~ 1 /im (HWHM) [28 . The 
laser (wavelength 532 nm) is detuned far off-resonance 
(wavelength 780 nm) , yielding a purely conservative and 
repulsive potential. The disorder amplitude Vr is ho- 
mogenous to better than 1% over the atom cloud (profile 
of half- widths 1.2 mm along y,z, 1 mm along x). 

The anisotropy of the speckle autocorrelation function 
(elongated along x) allows us to operate in a quasi-two 
dimensional configuration by launching the atoms per- 
pendicularly to the x-axis (along the z-axis). In the y-z 
plane, the atoms are scattered by a potential with a corre- 
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FIG. 3. Normalized angular profiles n(6,t)/n corresponding 
to the distributions shown in Fig. [2] The red solid lines cor- 
respond to the double structure fit around the backscattering 
direction [a parabola for the incoherent background (dashed 
line) and a gaussian for the CBS peak, see text]. Inset: false 
color representation of the momentum distribution (t = 2 ms) . 
The angle refers to the scattering direction with respect to 
the initial direction (pi). 



lation length shorter than the reduced atomic de Broglie 
wavelength (ki(j± c± 0.9), so that the scattering proba- 
bility is quasi isotropic. The dynamics within this plane 
develops on the typical time scale of a single scattering 
event, that is, the elastic scattering time r s . In contrast, 
the correlation length along the x axis is larger than the 
reduced atomic de Broglie wavelength (k\cr x ~ 4.5), so 
that each scattering event produces but a small devia- 
tion out of the y-z plane [29 . The diffusive dynamics 
along x is then slower than in the y — z plane, and at 
short time the dynamics is quasi 2-d. 

Figure [2] shows the time evolution of the momentum 
distribution for a disorder amplitude Vr/Ji = 975 ±80 Hz. 
In order to analyze these data quantitatively, we perform 
a radial integration of the 2-d momentum distribution on 
a thin stripe between p x — Ap res and p x + Ap res (inset of 
Fig. [3| (Ap res = mAv res is the momentum resolution). 
This yields the angular profile n (#,£), displayed in Fig. [3] 
for increasing diffusion times. 

We first extract the elastic scattering time r s from 
the exponential decay of the initial peak n(0 = 0,£) ex 
e~^ Ts [30]. We find r s = 0.33 ± 0.02 ms (mean free path 

= ^i r s = 1-1 /im). Note that this quantity also plays a 
role in the radial width of the ring (inset of Fig. |3| which 
is associated to a Lorentzian energy spread AE^is = ^/ r s 
(HWHM) acquired by the atoms when the disorder is 
suddenly switched on. Combining the corresponding mo- 
mentum spread with the resolution of our measurement, 
we find a width in agreement with the observed ring 
width. The measured value of r s is in accordance with 



numerical simulations adapted to our configuration, but 
is about 2 times the value predicted by a perturbative 
calculation [28]. This is consistent with the fact that 
we are not fully in the weak disorder regime defined by 
AE dis /E K = 2/k{l 8 <C 1. Here we have AE dis /E K ~ 0.4 

(Ms -5) [31]. 

Monitoring the isotropization of the momentum dis- 
tribution, we obtain another important quantity: the 
transport time r*, after which, in the absence of coher- 
ence, the information about the initial direction would be 
lost. It is determined from the exponential decay e~ l l T 
of the first component of the Fourier series expansion of 
n(0,i) [20] [33]. We find r* = 0.4 ± 0.05ms, also in agree- 
ment with numerics. Note that r* takes into account the 
CBS phenomenon and its calculation must include weak 
localization corrections. It is only slightly larger than r s , 
as expected for a nearly isotropic scattering probability 
(ha± = 0.9). The transport time sets the time scale of 
the onset of the diffusive dynamics, which is well estab- 
lished only after several r*. In the experiment, we ob- 
serve that the momentum distribution has become fully 
isotropic after t ~ 2.5 ms (i.e. ~ 6r*), with a steady 
and fiat angular profile of mean value n(#, t) = n, except 
around 6 = tt where the CBS peak is still present. 

To analyze the evolution of the CBS signal, we 
fit (see Fig. [3| the angular profiles by the function 
^incohOM) + n coh (£)exp[-(# - 7r) 2 /2A#(f) 2 ], where we 
assume that ni nco h(#,£) has a parabola shape around 
= 7r. This allows us to measure the contrast C(t) = 
^coh(£)/^incoh(7r, t) and the width A0(t) of the CBS peak. 
Their evolutions are shown in Fig. [4] A CBS peak ap- 
pears as soon as scattering in the backward direction 
is significant, but the contrast starts decreasing before 
reaching the ideal value of 1. For the width, we observe 
the predicted monotonous decrease, but it tends asymp- 
totically towards a non null value rather than zero. 

To compare these observations with theoretical pre- 
dictions, we must take into account our finite resolution. 
The black line in figure [4^b) represents the calculated 
CBS width that results from the convolution of our res- 
olution A# res = Ap res /pi = 0.07 with the expected CBS 
width A#cbs = ApcBS,a/Pi = h/piV2Dt in the diffusive 
regime (the widths are added quadratically) . The dif- 
fusion constant is evaluated using the standard relation 
D = vfr* /2, so that the solid line does not involve any 
adjustable parameter. We see that the agreement with 
the data is good when we enter the multiple scattering 
regime (for t > 1.5 ms in Fig. Qb)), but not at short 
times. 

The broadening of the CBS peak by the finite resolu- 
tion is also responsible for a decrease of the contrast, as 
represented in Fig. [4^a). Here also, we observe that the 
prediction (which again involves no adjustable parame- 
ter) is very different from the observed values at short 
times, but is in reasonable agreement with the measure- 
ments around t ~ 1.5 ms. On the other hand, the mea- 
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FIG. 4. Time resolved dynamics of the CBS peak: evolution 
of the contrast C (a) and the width AO (b) versus the scat- 
tering time t. The blue points are experimental datas, the 
solid and dashed black lines are the theoretical predictions 
for the multiple scattering regime at long times (t > 4t*; 
t* = 0.4 ms). The dashed and dotted red lines corresponds 
to the calculation of [34] at short times (i.e. for t ~ r s ; 
t s = 0.33 ms), where single scattering events cannot be ne- 
glected. 



sured contrast is definitely smaller than the theoretical 
prediction when t increases yet more. We relate this ob- 
servation at long times to the onset of the dynamics in 
the third direction x. Using a second imaging system 
yielding the momentum distribution in the x — y plane, 
we estimate a typical time of 4 ms for this out-of-plane 
dynamics to become significant, and render the 2-d ap- 
proach wanting. 

Deviations at short time were to be expected. In- 
deed, CBS demands multiple scattering, or at least dou- 
ble scattering, to happen (see inset of Fig. [I]), whereas 
single scattering events do not participate to the CBS 
peak. At short times (t ~ r s ), the contribution of sin- 
gle scattering to backscattering is not negligible com- 
pared to multiple scattering. This entails a reduction 
of the contrast (see e.g. [3 ), and a modification of 
the shape (no longer Gaussian), whose width decreases 
at this stage as 1/t (ballistic motion between the first 
two scatterers). In the case of light, a calculation for 
isotropic scattering [34] predicts a short time evolution 
of the contrast C = (2t/nr s )/{l + 2t/7rr s ) and width 
A#cbs ^ 3/kil 8 (r 8 /t). This prediction is plotted in Fig. [4] 
and is found in fair agreement with the observations in 
this time domain. Finally, note that the width around 
t ~ r s is linked to disorder strength quantified by k{l s . 
Here we find a maximum value of A# max ~ 0.3 rad, that 
is A# max ~ l.b/k[l s (k[l s ~ 5, see above). 

Similar measurements and analysis have been repeated 
for weaker disorder (Vr/H = 525,750 Hz) and smaller 
initial momentum pi {E^/h = 160,220,620 Hz), and we 
have found a similar agreement between data and the- 



ory. In contrast to the observed moderate changes in the 
maximum peak contrast (C max ~ 0.5— 0.7), the maximum 
peak width A# max increases significantly with the ampli- 
tude of the disorder and the inverse of p\. The highest 
observed width of 0.7 rad (from which we infer k[l s ~ 2) 
suggests that we are close to the strong disorder regime, 
where AL is expected to be experimentally observable in 
2-d systems. Such an observation, however, would de- 
mand a longer 2-d evolution in the disorder, which is 
limited in the present experiment because of the cross 
over to the 3-d regime. Increasing this time, as for in- 
stance in [35] , will then constitute the next step towards 
AL, with the possibility to observe the coherent forward 
scattering peak predicted in [36] . 

In conclusion, we have demonstrated experimentally 
that the time resolved study of the momentum distri- 
bution of ultra-cold atoms in a random potential is a 
powerful tool to study quantum transport properties in 
disordered media. We have been able to extract the elas- 
tic scattering time, the transport time, and to observe 
the CBS peak. Let us emphasize that the theoretical 
analysis as well as numerical simulations render an ac- 
count of the observations not only in the multiple scat- 
tering regime but also at short time, during the onset 
of multiple scattering. Such agreement gives a strong 
evidence of the fundamental role of coherence in that 
phenomenon. Further evidences of the role of coherence 
could be seeked in the predicted suppression of the CBS 
peak [37], when scrambling the disorder, or when dephas- 
ing the counter-propagating multiple scattering paths us- 
ing artificial gauge fields [38], in the spirit of pioneering 
works in condensed matter physics [13] or optics [39] . Fi- 
nally, this work also opens promising prospects to study 
the effect of interactions on CBS (see e.g. f4Q| [4Tj ) . 
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